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ABSTRACT 


A study of the effectiveness of atomic oxygen and 
carbon monoxide as cooling agents in a planetary thermo- 
sphere is carried out. Analytic expressions for the 
total flux radiated by these constituents are derived 
based on the assumption of a sharply defined radiation 
zone in which the scale height variation is linear. It 
is found that the total flux radiated by atomic oxygen 
can exceed that radiated by carbon monoxide if the 
temperature in the radiating region is less than 300°K. 

The time-dependent thermal conduction equation is 
considered briefly. It is shown that the cooling time 
for carbon monoxide is of the order of a few hours, in 
contrast to the several-day cooling time of atomic 
oxygen. A substantial overnight cooling is thus to be 
expected in planetary thermospheres which contain carbon 


monoxide . 
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1. Introduction 

Atomic oxygen and carbon monoxide are present in 
varying degree in the upper atmospheres of the Earth, 

Mars and Venus. The Earth's thermosphere is predomi- 
nately atomic oxygen and does not contain much carbon 
monoxide. The atmospheres of Mars and Venus, on the other 
hand, contain substantial CO^ and must therefore contain 
an abundance of both of the dissociation products, CO 
and 0 . 

Each of these constituents can cool a planetary 

thermosphere by infrared emission of some of the thermal 

energy deposited by solar extreme ultraviolet (EUV) 

radiation. The possible importance of the magnetic 

3 3 

dipole transition 0( P^) -* 0( P^) within the ground state 
triplet of oxygen in cooling the Earth's thermosphere was 
pointed out by Bates (1951). McElroy et al. (1965) have 
studied the effectiveness of CO rotational transitions in 
cooling the Martian thermosphere. 

The power radiated by atomic oxygen is only weakly 
temperature dependent, whereas the power radiated by 
carbon monoxide is proportional to the square of the 
temperature. However, direct comparison of the power 



radiated by each constituent can be misleading. The 

physically significant quantity is the total energy loss 
2 

per cm column of atmosphere by each constituent. The 
calculations presented below show that the total energy 
loss by O and CO is more strongly temperature dependent 
than is the power radiated locally. The importance of 
the 62 u radiation of atomic oxygen relative to rotational 
loss by CO is significantly greater than could be infer- 
red solely from comparison of the power radiated by each 
constituent. If the temperature at the altitudes where 
radiative loss occurs is less than 300°K, atomic oxygen 
can be the dominant radiator. A precise determination 
depends on the scale height gradients in the radiation 
zone, the temperatures at which radiative loss commences 
(assuming it commences at a well-defined optical depth 
in each constituent) , and whether or not the constituents 
are in diffusive equilibrium. Even for temperatures in 
the radiation zone of ~500°K the ratio of 0/C0 energy 
loss may be ~.5. 

The analytic expressions developed in this paper for 
the flux radiated by O and CO may serve as a useful check 


on the results of detailed integrations of the thermal 
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conduction equation. 

The time-dependent form of the thermal conduction 
equation is considered with a view to obtaining a quali- 
tative idea of the amplitude of the temperature drop 
expected during the night. When the terms representing 
conductive loss and radiative loss due to atomic oxygen 
are dropped, a simple expression for the temperature 
decay may be obtained. A cooling time for CO radiation 
is defined as the time in which the temperature would 
fall to half its initial value due to CO radiative loss 
alone. This cooling time is found to vary from about 
7 hours for an initial temperature of 1000°K to 12 hours 
for an initial temperature of 600°K. In the absence of 
a nocturnal heat source, CO cooling alone will cause the 
exospheric temperature to fall to £300°K in a 12 hour 
period for initial temperatures ;C1000°K. 

2. Power Radiated Locally by O and CO 

Radiative loss via transitions within the ground 

state triplet of atomic oxygen were considered by Bates 

(1951) in conjunction with his discussion of heating of 

the Earth's upper atmosphere. The ground state of 

3 3 

atomic oxygen is an inverted triplet with the Pq and P^ 
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levels lying .028 and .020 e.v. respectively above the 

P ^ level. The Einstein coefficients for these transi- 

-5 -1 -5 -1 

tions are A Q2 = 1.7 x 10 sec and A^ 2 = 8.9 x 10 sec 

Due to the greater excitation energy of the level, 

the small statistical weight, g Q = 1, and the relative 

magnitudes of the Einstein coefficients, emission from 

3 3 

Pq will be much smaller than that from P^ and may be 

neglected. The question of thermodynamic equilibrium 

within the 0(^P ) triplet was considered by Bates, 

\J , L t Z 

who concluded that collisional excitation of the excited 

levels would probably occur much faster than photon 

9 -3 

emission if the total particle density were ^.10 cm 
It will be shown below that infrared emission in the 62u 
oxygen line is almost certain to be effective at much 
higher densities and that the assumption of equilibrium 
within the oxygen ground state is therefore justified. 

With this assumption the power radiated in the 62u line 
is 
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3 3 

where ^ and e are the energies of the and levels, 

g^, g^, and g^ are the statistical weights (the subscript 
refers to the corresponding J value) and n refers to the 
oxygen number densities. Evaluating the constants results 
in 


R 

s 


2.85 x 10 


-18 


n ox (z> 


f ( T) ergs/cm sec 


(2) 


where is the total oxygen density and 

CJX ! 


f (T) = 


. -228/T 
3e 
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The temperature dependence of the power radiated locally 
by atomic oxygen is not very great. The factor f(T) varies 
from 0 to 1/3 as T varies from 0 to 00 . 

The expression for the power radiated locally by CO 
was also given by Bates. The power radiated is proportional 
to the square of the temperature in this case and is given 
by 


R 

r 


-23 


n 


CO 


2.58 x 10 


T^ ergs/cm^ sec 


(4) 
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Rotational states are readily populated by collisions 
(Zener, 1931) and it is highly probable that carbon 
monoxide is in rotational equilibrium well above the 
region of radiative loss. 

3. Equation of Radiative Transfer for a Doppler -Broadened Line 
An accurate evaluation of the total energy lost by 
a radiating constituent such as 0 or CO requires a deter- 
mination of the altitude at which unit optical depth 
occurs. We will assume that radiative loss is fully 
effective above unit optical depth and ineffective at 
lower altitudes. An expression for the optical depth in 
Doppler -broadened lines may be obtained from the equation 
of radiative transfer. We consider a two-level system 
with the upper energy state denoted by subscript 2 and 
the lower state by subscript 1, and a radiation field of 
density f.: . The rate of decrease in the population of 

the upper energy state is 


dn 2 

at- = n 2 A 21 + n 2 B 21 C v - n l B 12 e v 


( 5 ) 
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spontaneous and induced emission and for absorption. 

Each transition resulting in a decrease in corresponds 

to the addition of a photon to the radiation field. Since 

this emission occurs over a frequency range the decrease 

in n 2 corresponds to an increase in the photon intensity 
-2 -1 -1 -1 

(in units cm steradian sec (cps) ) integrated over 
frequency and direction, i.e., 


dn. 


- at 2 ■ I dn a§J' ' ^ t 


( 6 ) 


where dQ is the element of solid angle, I the radiation 
intensity, and ds an element of path length in the medium. 
If the spectral line due to transitions between states 
2 and 1 is Doppler-broadened we may write 
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where a n is the Doppler width and I the intensity at the 

V 0 

line center. We now make use of the relationship between 
the radiation density and intensity, 
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between the Einstein coefficients 
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21 

L 21 
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and 
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where and are the statistical weights of levels 2 
and 1, and substitute the combined results of (6) and (7) 
in (5) to obtain the radiative transfer equation 
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where F., = - 

V TT' 


yl do is the photon flux. The subscript 0 


on v has been dropped. The relation between the direc- 
tions specified by z and s is dz = ds cos 0 where 9 is 
the angle between directions z and s. If we now assume 
that the levels are in thermal equilibrium we have 


n _l = g _l e 
n o 9o 6 


, 21 /kT 


and 


n 2 Z(T) 


2 - £ 2 AT 
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where = ^2 ~ e l the energ y difference between 

states 2 and 1, n is the total number density and Z(T) 
is the partition function for the system. With the use 
of (11) the radiative transfer equation is 
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where 
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is the absorption cross section in the line and 
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is the emission coefficient. The subscript D refers to 
the fact that these expressions hold for Doppler-broadened 
lines. The ratio e D / a D is just the Plank distribution 
(for photons) as it must be since thermal equilibrium has 
been assumed. 

4. Optical Depths in 0 and CO 

We will obtain some initial estimates of the over- 
lying number densities in oxygen and carbon monoxide at 
which these constituents may begin to cool a planetary 
thermosphere by radiative loss. The optical depth in a 
line is defined by 


d'^ = — n n dz 


(15) 
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where the altitude variable is denoted by z . If the con- 
stituent n is distributed in diffusive equilibrium in an 
isothermal atmosphere, the optical depth is 


t = a D n H 


( 16 ) 


and this unit optical depth occurs at an altitude where 
the overlying number density is 


n H = — (17) 

°D 

For the 62u line we evaluate the constants in (13) and 

3 

use the partition function appropriate to the 0( P ) 

U t X / £* 

triplet to find 


ct d ( 0) = 8.9 x 10 
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If we take T = 300°K as an estimate of the isothermal 
temperature, then the overlying oxygen number density at 
unit optical depth is 


n 


H 

OX OX 1 T=1 


= 2.4 x 
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We recall that the initial assumption of thermal 

equilibrium in the ground state oxygen triplet required 

9 

total number densities ^10 . That radiative cooling 

commences below the altitude at which departures from 

equilibrium occur is obvioas from (19) since atomic 

9 

oxygen number densities of 10 would correspond to unit 
optical depth in the 62u line only if the scale height 
were -1000 km. 

In the case of carbon monoxide we have transitions 
between rotational states J ^ J-l for all values of the 
angular momentum quantum number J. We will assume, as 
has been done by McElroy et al. (1965), that the optical 
depth in the carbon monoxide rotational lines is given 
by the optical depth in the line in which most energy is 
transferred. The rotational state with the greatest 
occupation is specified by 


J 


= r T 

MAX ! 2 $~ 


1/2 


( 20 ) 


The Einstein coefficient for spontaneous emission in this 
line may be evaluated from the general expression for the 
A coefficient (Schiff, 1955, Chapter 10) using the dipole 
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matrix element for adjacent states of a rigid rotator. 
We have 


J, J-l 


4 , Aw ^,3^2 J 
3h (4 r ) d 2J+1 


( 21 ) 


where B is the rotational constant of the molecule, d 
the dipole moment of CO and R Plank's constant. Evaluat- 
ing this and all other J-dependent terms in (13) at 
J = gives the absorption cross section in the 

strongest rotational line as 


D 


(CO) 
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where is the characteristic rotational temperature 

— 18 

hcB/k. Evaluating the constants, with d = . 12 x 10 - 
(Fowler, 1935), B = 1.93 cm 4 (Bates , 195 1) , and expanding 
the exponential to first order we find 


.(CO) = 


1.35 x 10 

ZJ7T 


-13 


cm 


( 23 ) 


This expression was given by McElroy e_t al. (1965). At 
a temperature T = 300°K unit optical depth in the CO 
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rotational lines thus occurs at an overlying CO density 
of 

n H i . = — = 3.8 x 10 16 cm" 2 (24) 

CO CO 1 t=1 (CO) 

Equations (19) and (24) show clearly that conclusions as 
to the relative effectiveness of 0 and CO as cooling 
agents based on a direct comparison of the power radiated 
locally is misleading. The sum of the oscillator strengths 
for the CO rotational lines is greater than the oscillator 
strength of the 62;j 0 emission. However, this implies 
not only that CO is a more effective emitter in the infra- 
red than is O, but is also a more effective absorber. 

Only a consideration of the balance between emission and 
loss can lead to valid conclusions regarding the relative 
effectiveness of these constituents as radiators. 

5. Comparison of Total Fluxes Radiated by O and CO 

We will now assume that in the altitude region 
from which radiative loss occurs the scale height may 
be adequately described in terms of a linear variation 
i . e . , 
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H(z) = H + 0 z 


(25) 


where is the scale height at the bottom of the zone 
and b is the (constant) scale height gradient. The 
number density of a constituent is given in this case by 


n (z) = 


n T 
0 0 

tTzT 


1/8 


*‘0 

H 


(26) 


The subscript 0 on a quantity refers to the value of 
that quantity at the altitude of the bottom of the radia- 
tion zone, i.e., at unit optical depth. This altitude 
is different for different constituents, but this dis- 
tinction need not be made yet. The optical depth in a 
given constituent is 


(z o } = 


n(z) 




dz 


'0 


(27) 


where we assume that z^ denotes the altitude of the top 
of the radiating zone and that contributions to the opti- 
cal depth from greater altitudes are negligible. The 
altitude z^ is not strictly defined, but should be 
chosen at least one scale height above Zq. Using (18) 
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and (23) for the absorption cross sections in 0 and CO 
we find the following expressions for the optical depths 
in these constituents 
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(29) 


where a = T^/T^ is the ratio of temperatures at the 

bottom and top of the radiation zone. Subscripts OX and 

CO have been appended to quantities which refer to 

atomic oxygen and carbon monoxide. In the derivation of 

— 1/2 

(28) from (27) only the T dependence of - (T) was 

explicitly included in the integration. The remaining 
temperature-dependent term in (18), which varies from 0 
to 1/5 as T varies from » to 0, is assumed evaluated at 
Tq . This approximation will slightly overestimate the 
optical depth at 62u (for a given altitude) and corre- 
spondingly underestimate the effectiveness of O cooling. 
The total flux radiated by a constituent is given by 
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R(z,T) 


dz 


(30) 
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where R(z,T) is the power radiated. Using (2) and (4) 
for the power radiated by 0 and CO we find 


f r(°) 

f r (co) 


— 18 r 

= 2 * 8 x 10 ^oVox^Vt 1 - 
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The function f(Tg) appearing in (31) and given explicitly 
by (3) was not included in the integrand of (30) but is 
assumed evaluated at T^. The variation of f(T) with T 
is relatively slow for T > 250°K. The expression (32) 
was given by McElroy et al^ (1965). If we now evaluate 
the product nH appearing in (31) and (32) at unit optical 
depth (using (28) and (29)) we obtain the total energy 

loss by O and CO in a planetary thermosphere (in ergs/ 

2 s 

cm sec) 
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where 


E( V 228/T q 

e - 1 

It must be borne in mind that T Q may refer to different 

temperatures in the expressions relating to 0 and CO. 

For temperatures T Q 300°K the flux lost in the 62u 

3/2 

line is roughly proportional to T Q since E(T q ) T Q . 

The presence of atomic oxygen in a planetary thermosphere 

therefore does, to some extent, act as a thermostat since 

the total energy loss increases more rapidly than the 

temperature. In an atmosphere containing both 0 and CO 

the thermostatic effect of CO will of course predominate 

7/2 

since the energy loss by CO goes as 

We will neglect the difference in tempera lures at 
unit optical depths in 0 and CO and denote by T R a mean 
temperature in the altitude region between which 0 and 
CO begin to radiate. The ratio of the total heat 
fluxes radiated by 0 and CO is then 
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Since E(T ) ~ T for T 300°K this ratio falls off 
R R R 

-2 . 
roughly as T at high temperatures. The ratio is 

plotted as a function of T in Fig. 1. Table 1 gives 

the values of the total flux radiated by O and CO from 

unit optical depth in each constituent for various 

values of temperature and scale height gradient in the 

radiation zone. The figures may serve as a useful check 

on the results of detailed numerical integration of the 

thermal conduction equation and as a basis for judging 

the relative effectiveness of these constituents as 

cooling agents. The temperature difference at unit 

optical depth in O and CO should be taken into account 

before quantitative comparisons are made. 

6. Nocturnal Temperature Change 

Bates compared the total energy content of the Earth's 
2 

F region (per cm column) with the energy loss rate due 
to thermal conduction and 62u O radiation. The character- 


Total flux radiated from a planetary thermosphere by O and CO. T refers to the 
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istic time which he obtained as the quotient of these 
quantities was of the order of several days, and he 
therefore concluded that nocturnal cooling would be 
slight. The time-dependent study by Harris and Priester 
(1962) gives a daily amplitude of exospheric temperature 
variation of roughly ^200°K about a mean of 1500°K. The 
amplitude is thus relatively small (-15-20%) . 

The situation is quite different if substantial CO 
is present. We will define the cooling time for CO as 
the time required for radiation from the rotational 
states to halve a given initial temperature. If we con- 
sider CO cooling alone, the temperature change is 
described by 


dT 
pC p dt 


n co (z) T 


(37) 


23 o 2 

where k = 2.58 x 10 erg/ K sec is the coefficient 
r 

of the term for the power radiated by CO (see equation (4)) . 
The solution of (37) is 
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where T Q is the initial temperature and 
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(39) 


is the cooling time for CO. If we take the density p in 
(39) to be the sum of the 0 and CO densities, 0 may be 
expressed as 


, C P m COr , 
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where m__ is the mass of a CO molecule, m the mass of 

CO ox 

atomic oxygen, - D (0) and o D (CO) are the Doppler absorption 
cross sections given in (18) and (23) , and and r CQ 

are the isothermal optical depths in O and CO. The 
temperature dependence not taken into explicit account 
in passing from (37) to (38), as well as the altitude 
dependence of 8 , appears in the right hand term in the 
brackets of (40) . The Doppler absorption cross sections 
are temperature-dependent, but the ratio ct^(C0)/g^(0) 
does not depend strongly on temperature and is equal to 
about 5 over a wide temperature range. The first two 
factors in the right hand term in brackets of (40) are 
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thus of order unity. The ratio 'qx^ CO '*" s certa ^ n -*-y less 
than one, but since the density, V , is likely to be 
slightly underestimated in (40) we will take the right 
hand term as equal to 1. The cooling time is thus 
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c m__ 
p CO 

k T 
r 0 


(41) 


The approximation made in setting • q%/ ' 1 is thus 

roughly equivalent to the assumption that the carbon 

3 

monoxide density (gm/cm ) in the radiation zone is half 
the total density. Equation (41) thus contains the 
implicit assumption that CO^ is a major atmospheric 
constituent at the altitude at which it dissociates. 

The temperature change as a function of time, eval- 
uated from (38), over a 12 -hour period is shown in Fig. 2 
for various initial temperatures. The specific heat 
has been taken as c^ = 0.2 cal/°K gm. The cooling time 
given by (41) is about 7 hours at 1000°K. 

The amplitude of the temperature decrease given in 
Fig. 2 should be representative of that to be expected 
in the Martian thermosphere since the rotational period 
of the planet is about 24 hours. Thus, if the peak 



24 


daytime temperature is £1000°K the temperature of the 

o 

nighttime Martian exosphere must fall to ^300 K unless 
there is a nocturnal heat source. The rotational period 
of Venus is relatively long, being about 250 days retro- 
grade (Goldstein and Carpenter, 1963). The nocturnal 
temperature in the Venusian thermosphere could thus be 
quite low since equilibrium may depend either on energy 
transport from the day side, on infrared heating from 
the planetary spectrum, or on other low altitude heat 
sources. Thermal conduction from the day side could 
not maintain a very high temperature on the night side 
of the planet. Even with a day-night temperature dif- 
ference of 2000°K, the mean horizontal temperature gradient 
would be . 05°K/km, which would transport ~.008 erg/cm^sec. 
From Table 1 we see that this flux could be substantially 
accommodated by O and CO at temperatures --■ 100°K. A 
quantitative estimate of the temperature maintained by 
infrared heating would require a solution of the radiative 

transfer equation, but it would be less than the plane- 

o 

tary effective temperature of ~230 K. 

7. Summary and Conclusions 

The effectiveness of atomic oxygen and carbon monoxide 



as radiative cooling agents in a planetary thermosphere 
has been investigated. Analytic expressions for the 
total flux radiated by each constituent were developed 
and it was shown that for temperatures -•t300°K in the 
radiation zone the total energy loss in the 62 L ; oxygen 
line may equal or exceed that in the CO rotational lines. 
The figures given in Table 1 for the total flux radiated 
by each constituent depend on the efficacy of a Bates' 
type of cooling approximation for their validity, but 
the conclusions drawn in this paper as to the relative 
cooling effectiveness of O and CO are not so dependent. 

It has been shown that substantial overnight cooling 
will occur in planetary exospheres containing CO. The 
amplitude of the day-night temperature variation in the 
Martian exosphere is probably about 50% in contrast to a 
15-20% diurnal variation in the Earth's exosphere. The 
nocturnal temperature minimum in the Martian exosphere 
must be <300°K if the peak daytime temperature is ^1000°K, 
as it must be for the retention of atomic oxygen by the 
planet (Chamberlain, 1962) . The Venusian exosphere must 
also be rapidly cooled by CO. Due to the length of the 
night on Venus the equilibrium exospheric temperature 


must be maintained either by horizontal energy transport 
or by heating from lower altitudes. Horizontal conduc- 
tion from the day side could not maintain a nighttime 
temperature in excess of — 100°K. 
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LEGENDS 


Fig. 1. Ratio of 0 to CO radiated flux. The 
curves are labeled according to the values of the scale 
height gradient fl for 0. 

Fig. 2. Decrease in temperature with time due to 
infrared radiation from CO. 
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